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The appearance of the young supernova remnant SN1006 is dominated by 
emission from non-radiative shocks in the NE and NW regions. At X-ray ener- 
gies the northeast shock exhibits predominantly nonthermal synchrotron emis- 
l/^ ' sion, while the northwest shock exhibits a thermal spectrum. We present far 

ultraviolet spectra of the northeast (NE) and northwest (NW) portions of SN 
1006 acquired with the Far Ultraviolet Spectroscopic Explorer (FUSE). We have 
detected emission lines of O VI (1032, 1038 A) and broad Ly-fS (1025 A) in the 
^ ■ NW filament, but detect no emission lines in the NE region down to a level of 4.7 

■ X 10^^'' erg cm~^ s~^ arcsecond"^. We observed in the NW an O VI intensity of 

^ ' 2.0 ± 0.2 X 10^^^ erg cm~^ s~^ arcsecond"^ and measured an O VI line width of 

O ■ 2100 ± 200 km s~^ at a position where the Ha width was measured to be 2290 

^ ■ ± 80 km s~^ (Ghavamian et al. 2002). This implies less than mass proportional 

heating of the ions. Using the ratio of intensities, I(NW)/I(NE)~n(NW)/n(NE), 



2 ■ the density ratio of the two regions is found to be > 4, a value that is with the 

■ uncertainties of the ratio of 2.5 measured by Long et al. (2003). The derived 

O VI kinetic temperature is compared to previous estimates of electron, proton. 



and ion temperatures in the remnant to study the relative heating efficiency of 
various species at the shock front. The degree of postshock temperature equi- 
libration may be crucial to particle acceleration since the temperature of each 
species determines the number of high speed particles available for injection into 
an acceleration process that could produce Galactic cosmic rays. 

Subject headings: ISM:individual(SN1006)-supernova remnants-shock waves- 
ultraviolet:ISM 
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1. Introduction 

SN1006 (G327.6+14.6) is a nearby Type la supernova remnant at a distance of 2.1 kpc 
(Winkler, Gupta, Long 2003). With a mean expansion rate of 8700 km s~^ it is ~18 pc 
wide (Winkler, Gupta, Long 2003). The remnant has a high Galactic latitude and modest 
foreground reddening, E(B-V)=0.11 ± 0.02 (Schweizer & Middleditch 1980). This young 
supernova remnant is entering the Sedov- Taylor phase of supernova remnant evolution. 

SN1006 has been observed at radio (Pye et al. 1981), optical (Ghavamian et al. 2002; 
Kirshner, Winkler, and Chevaher 1987; Smith et al. 1991), ultraviolet (Raymond, Blair, 
Long 1995) and X-ray (Winkler, Gupta, Long 2003; Long et al. 2003; Bamba et al. 2003) 
wavelengths. Gamma ray observations (Tanimori 1998) were reported but not confirmed. 
Thin, pure Balmer line filaments were found in the optical. In the radio and X-ray, the 
remnant has a limb-brightened shell structure with cylindrical symmetry around a SE to 
NW axis probably aligned with the ambient galactic magnetic field (Reynolds & Gilmore 
1986; Jones & Pye 1988). The NE shock front of SN1006 shows strong non-thermal X-ray 
and possible gamma ray emission while the NW shock shows very little non-thermal emission 
at radio or X-ray wavelengths. 

Ly-/9, He II, C IV, and O VI lines were observed from the faint optical Balmer line 
filament of the NW shock of the supernova remnant, by the Hopkins Ultraviolet Telescope 
(HUT), flown during the Astro-2 space shuttle mission. The observed FWHM of the lines 
were 2230, 2558, 2641 km s~^, respectively (the O VI line width could not be measured). A 
kinetic temperature could be calculated from these line widths. The kinetic temperatures of 
these species are not equal, because they do not scale inversely with the square root of their 
atomic mass. Instead, the UV observations do suggest that Tion ~ ^^^Tproton indicating 

Hip 

little to no temperature equilibration between species. 

SN1006 provides an opportunity to investigate parameters of non-radiative coUisionless 
shocks faster than 2000 km s~^. CoUisionless shocks appear in many astrophysical phenom- 
ena, from coronal mass ejections (CMEs) in the hehosphere to jets in Herbig-Haro objects. 
When a shock is non-radiative the detection of emission from the shock front is possible, 
as all of the optical and UV emission of a non-radiative shock comes from a narrow zone 
directly behind the shock front. Interactions at the coUisionless shock front depend upon 
mechanisms such as plasma waves to transfer heat, kinetic energy and momentum, and it 
is not well understood how particles of different masses and charges are affected by these 
processes. The temperature of the species and the degree of temperature equilibration be- 
tween electrons, protons and other ions is central to the interpretation of X-ray spectra, 
which effectively measure electron temperature. The energy distribution of a particle species 
is important to cosmic ray studies as only those particles at a high energy tail of a particle 
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distribution are available for cosmic ray acceleration. 

The method of using Ho; hnes to determine coUisionless shock parameters was originated 
by Chevalier & Raymond (1978) and Chevalier, Kirshner, & Raymond (1980). The Ho; hne 
has a two component profile. The width of the broad component of the Ha line is related 
to the post-shock proton temperature as a result of charge exchange between neutrals and 
protons, which produces a hot neutral population behind the shock. The narrow component 
of the Ha line is produced when cold ambient neutrals pass through the shock and emit line 
radiation before being ionized by a proton or electron. The ratio of the broad to narrow flux 
is sensitive to electron-ion equilibrium and the pre-shock neutral fraction. The FWHM of 
Ha line was measured to be 2290 ± 80 km s^^, with models implying the speed of the shock 
is Vghock = 2890 ± 100 km s"^ (Ghavamian et. al. 2002; hereafter GWRL02). The Ha broad 
to narrow intensity ratio measured to be 0.84 implies an electron temperature much lower 
than the ion temperature. 

This UV observation from the FUSE satellite focused on the shock front in the NW 
observed by HUT and GWRL02 and on a region in the NE dominated by non-thermal 
emission. Prom the spectra, a broad Lyman /3 line (A~1025 A) and the doublet of O VI 
(A~1032, 1038 A) were analyzed for spectral width, intensity, and flux. We use the line 
widths of the NW and the intensities of the O VI lines in the NE and NW shock fronts 
to compare the electron-ion and ion-ion temperature equilibration efficiencies as well as 
densities. The heating of different particle species by the shock front as well as parameters 
of coUisionless shocks that affect particle species heating will be discussed. 

2. Observations 

The Far Ultraviolet Spectroscopic Explorer (FUSE) has a wavelength range of approx- 
imately 900-1180 A. The Large Square Aperture (LWRS), with a field-of-view of 30" x 30", 
with a roll angle of 167°, was chosen for this observation because models predicted that the 
O VI emission behind the shock would be spread over 35" (Raymond, Blair, Long 1995; 
Laming et al. 1996). The LWRS has a fiUed-aperture resolution of about 100 km s^^. The 
FUSE spectrometer consists of four independent channels with two segments each. Four 
of these eight segments operate in the wavelength range for the O VI doublet, A=1031.91, 
1037.61 A. However, the Sihcon Carbon (SiC) coated channels, because they are optimized 
for A < 1020 A, add an unacceptable amount of noise to our faint signal, so only the Lithium 
Floride (LiF) channels are used. These two segments are designated LiFlA and LiF2B. The 
LiFlA channel covers wavelengths 987.1 - 1082.3 A, while the LiF2B covers 979.2-1075.0 A. 
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Although the northwest region of the remnant has been observed before in the UV 
(Raymond, Blair, Long 1995), we have much better spectral resolution and a more opti- 
mal aperture size to include the entire ionization region given that it may be larger than 
19"(Laming et al. 1996). The apertures used for past observations were HUT=19" x 197" 
(Raymond, Blair, Long 1995) and CTIO RC Spectrometer=2"x 51" (Winkler, Gupta, Long 
2003)(GWRL02). We positioned the aperture center to be 5"-10" behind the Ha filament 
where the peak formation of O VI occurs. The NE position was chosen based on the edge of 
the X-ray filament from Long et al. (2003). 

FUSE observations of the northwest region, centered at Q!20oo =15*^ 2™ 19. 17*^, ^2000 =- 
41° 44' 50.4" were obtained on 23 June 2001 and 26 February 2002 with total exposure times 
of 35,627 s and 6,690 s. Observations of the northeast region, centered at Q!2ooo=15^ 4°^ 5.0^, 
52ooo=-41° 50' 40.5" were obtained on 25 June 2001 and 27 February 2002 with exposure 
times of 42,365 s and 9,666 s. The locations of observations are shown superimposed on 
an Ha image of the remnant taken with the CTIO Schmidt telescope in Figure 1 (Winkler, 
Gupta, Long 2003). Inserted in the figure is a close up from Chandra (Long et al. 2003) of 
the NE region of observation to illustrate the x-ray morphology, although no optical emission 
is obviously present. 

There are four components to the background of this observation; detector background, 
geocoronal lines, the diffuse galactic UV continuum and diffuse galactic O VI emission. 
The background count distribution on the FUSE detectors is composed of two separate 
components (Anderson, Sankrit, Dupuis 2003). The 'intrinsic' background forms from the 
/3-decay of potassium in the microchannel plate (MCP) detector glass and the spacecraft 
radiation environment. The effect of the spacecraft radiation environment on the detector 
background varies from night to day and with solar activity, but over a short observing time 
this variation is not significant. The second component is caused by scattered light, primarily 
geocoronal Ly-a. This line produces detector averaged count rates as small as 20% of the 
intrinsic background during the night and increasing to 1-3 times the intrinsic rate during 
the day. The other two components of the background, galactic UV emission and diffuse 
O VI emission, will be discussed later. 

The observations were calibrated with the CalFUSE Pipeline Version 2.2.1. Data from 
all exposures arc processed through the pipeline and then co-added following the FUSE 
Data Analysis Cookbook and The FUSE Observer's Guide. The data were selected to 
contain only the night observations. This greatly reduces the geocoronal background. The 
night-only exposure times were 32,287 s for the Northwest and 39,386 s for the Northeast. 
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3. Analysis and Results 

As mentioned above, the background consists of detector noise, geocoronal lines, diffuse 
galactic O VI and an astrophysical UV continuum. The first two sources were explained in 
the previous section, but the additional diffuse UV continuum must be treated separately. 
It does not originate from SN1006, as it is seen in both of the entirely different regions of 
the remnant; the NE shock and the NW shock. The diffuse background is attributed to 
hght from hot stars scattering on dust. The diffuse UV continuum is especially bright in this 
region of the sky according to models by Murthy & Henry (1995). A value of 8.4 x 10~^^ erg 
cm~^ s~^ A~^ was quoted by Raymond, Blair, Long (1995) while we are seeing approximately 
6.5 X 10~^^ erg cm~^ s~^ through an aperture one quarter the size of the HUT observation. 

In addition to the diffuse UV continuum, Shelton et al. (2001, 2002) and Otte et. al. 
(2003) have found a diffuse O VI background. The brightness of the O VI background is 4700 
± 2400 photons cm~^ s^^ sr^^(Otte, Dixon, & Sankrit preprint). The widths of the diffuse 
O VI lines fall between 10 and 160 km s^^. In the current NE spectrum diffuse O VI emission 
has a width of < 200 km s~^ and a brightness of 3500 photon cm~^ s~^ sr~^. We attribute 
the NE emission to the diffuse galactic O VI background. This enabled us to subtract the 
NE as a background from the NW data to further ehminate airglow lines, the diffuse UV 
emission and the galactic O VI background. The intensities of the airglow lines at 1042A 
and 1048A are quite similar in both the NE and NW, further allowing this subtraction. The 
raw spectra of the NW and the NE regions are shown in Figure 2, with airglow lines marked. 

For the NW region, a nonlinear chi-squared minimization routine was used to fit Gaus- 
sian line profiles to the spectra. The wavelengths considered for analysis were restricted to 
1010-1050 A to minimize spurious background effects near the ends of the detector's spectral 
range. The data were binned by 0.1 A to increase the number of counts per bin without 
losing resolution, as the line widths were several Angstroms wide. The width of the broad 
Ha line, from GWRL02, is Vh=2290 km s~^. Since the Lj-f3 line is formed by the same 
process(Chevalier, Kirshner, & Raymond 1980), its line width was set equal to the Ha broad 
component width . The shift of the centroid of the broad and narrow component of Ha, 
v=29 km s~^, is effectively negligible (implying that the shock is viewed completely edge-on) 
so the broad Ly-/3 fine centroid was fixed at its rest wavelength. Only the intensity of the 
fine was a free parameter. The blue wing of the line was fit from 1010 A to 1024.5 A. Due 
to the extinction from interstellar dust, a correction factor must be applied to deredden the 
observed flux. Using the extinction curves of Cardelli, Clayton, Mathis (1989), the resulting 
dereddened Ly-/? flux is 2.3 ± 0.3 x 10~^^ erg cm~^ s~^. 

After subtracting the fitted broad Ly-/3 line profile, the wavelength range from 1022-1028 
A was excluded from the fitting routine in order to avoid negative fiuxes and residual airglow 
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that would skew the gaussian fits of the O VI hnes. At ~1037.0 A there were absorption 
features present that coincided with a C II hne and molecular hydrogen lines, along with an 
O I airglow line. The absorption feature with the spectral range from 1035 - 1038 A was 
therefore excluded from the fit. 

The O VI doublet was fit with two gaussians with fixed centers at 1031.91 and 1037.61 

A respectively corresponding to the ccntroid of Ha. The doublet line intensities were forced 
to have a 2:1 ratio but the magnitude of the intensities were allowed to vary. The observed 
flux is 6.7 ± 0.1 X 10^^^ erg cm^^ s~^ arcsec^^. Total dercddcncd flux for the O VI doublet 
lines was 1.8 ± 0.2 x 10^^^ erg cm~^ s~^. The O VI line widths were measured to be 7.2 ± 
0.4 A FWHM, or equivalently 2100 ± 100 km s^^. The formal error on the fit is 100 km s^^. 
However, due to systematic error a more conservative error of ± 200 km s"^ is used. The fits 
are shown in Figure 3. The width is within the limiting estimate of Raymond, Blair, Long 
(1995) of < 3100 km s"^ and is within la error of the Ha width of 2290 km s~^ Although 
the faint signal in the NE did not allow for a statistically significant fit, an upper limit of 
O VI intensity was found assuming a width of 2000 km s^^. The observed upper limit on 
the O VI line is 1.6 x 10^^^ erg cm^^ s^^ arcsecond^^. An upper limit on the dereddened 
intensity of O VI in the NE is 4.2 x 10^^^ erg cm~^ s~^. The upper limit of fiux for Ly-/? 
in the NE is 1.6 x 10~^^ erg cm~^ s^^ arcsecond"^. An upper limit to the dereddened Ly-/3 
intensity in the NE region is 4.4 x 10~^^ erg cm~^ s~^. 

Past observations of SN1006 line widths and intensities arc summarized in Table 1. In 
order to compare past measurements made with varying aperture sizes, we use intensity 
per arcsecond measured along the length of the filament. The Ly-/3 from the HUT and the 
current FUSE observation are consistent. We can use the various measurements to study 
the ion heating. The proton temperature was found using the shock speed of 2890 km s~^ 
from GWRL02. This proton temperature was then multiplied by mjo„/mp to calculate the 
mass proportional temperatures. These calculated temperatures were then compared to the 
temperatures given by using the FWHM of each ion line. The temperature of O VI as 
indicated by its FWHM is less than mass proportional by 48%. For the other ions. He II, 
C IV the heating was also less than mass proportional, by 21% and 18% respectively. 

The brightness of the O VI lines is proportional to density, uq, and the depth of the 
filament along the line of sight. Therefore, an upper limit to the density in the NE can 
be found by the ratio of intensities provided that the depths along the line of sight are 
known. Long et al. (2003) calculated a density ratio of n(NW)/n(NE) = 2.5. From the 
thermal component of the Chandra X-ray spectra Long et al. (2003) estimated a pre-shock 
ISM density of 0.25 cm~^ in the NW. Using the limit to the O VI intensity ratio of the 
NW and NE a ratio of the densities is found to be n(NW)/n(NE) > 4, which is within the 
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uncertainties of the Long et al. calculations. Therefore, assuming a pre-shock density in 
the ISM of 0.25 cm~^ in the NW, the pre-shock NE density < 0.06 cm~^. This density 
calculation depends on the assumptions of similar depths along the line of sight in the NE 
and the NW and of similar numbers of O VI photons per atom passing through the shock. 
The amount of electron-ion equilibration in the NE would affect these assumptions. Greater 
electron-ion equilibration in the NE would increase the number of O VI photons per atom 
(Laming et al. 1996), so the limit on the density in the NE would be even smaller. We 
attribute the low upper limit on the O VI intensity in the NE to the low density medium 
into which the remnant is expanding. 

4. Discussion 

O VI lines were not conclusively observed in the faint non-radiative non-thermal NE 
shock indicating that the two distinct shock regions heat ions differently. Ion heating is 
important to cosmic ray acceleration and the overall energy distribution of the system. 
The ions have most of their kinetic energy in a broad distribution which is generally non- 
Maxwellian as the time to equilibrium via Coulomb collisions for ions and protons is 1.2 x 
10^ years. To understand the heating at the shock front, turbulence, line widths, methods 
of calculating heating, and the role of neutrals at the shock front will be discussed. 

4.1. Small Scale Turbulence 

Turbulence plays a role in the evolution of fast shocks in supernova remnants (Reynolds 
2004; Ellison & Reynolds 1991). Small scale turbulence spreads the line profile of an ion 
much like thermal broadening of a line profile. Since some of the shock energy must be 
used for bulk flow, we will examine turbulence with a velocity of 1500 km s~^ which is large 
enough to affect the spectra but not contain all the energy of the flow. Turbulence decays 
on a time scale proportional to the characteristic length of the turbulence divided by the 
velocity of the turbulence ~ £/v (Tennekes & Lumley 1972). The width of the Ha filament 
is at most 10^^ cm based on its 1" apparent width on the sky(GWRL02), making the time 
scale of the turbulence 10^ s ~ 3 years. Using this decay time and the post-shock speed of 
750 km s~^, one quarter of the shock speed, the post-shock region affected by turbulence 
would be 7.5 X 10^^ cm. The O VI filament with an observed width of 3 x 10^''' cm, assuming 
the 30" FUSE aperture is filled, is also too wide to be dominated by turbulence thus the 
turbulence that is present in the shock of SN1006 is short-lived and not a major source of 
line broadening. 
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4.2. Line Widths of O VI, UV lines and Ha 

The UV line profiles of the current observations can be compared with past observations 
of various ion species. The currently observed O VI line width in the NW shock is within 
la of the Ha line width previously measured by GWRL02. Vink et al. (2003) measured 
an O VII line width of 3.4 ± 0.5 eV, or approximately 1775 ± 261 km from a different 
northwest region. This line width is substantially narrower than those of other ion species 
measured thus far, although the region of observation for this measurement is different from 
the position of our observations. Along a 124" slit, Smith et al. (1991) found little variation 
in the Ha profiles, indicating that the oxygen temperature does not vary significantly along 
the length of the NW filament. This implies one of two processes. First, the line width 
could decline with ionization state and distance behind the shock due to Coulomb collisions, 
as Coulomb collisions would transfer heat to other species. In Section 4.1, we found the 
Coulomb coUision time to be far too long for this process to be important. The second more 
probable scenario is that some of the lower temperature O VII is from the reverse shock in 
the supernova ejecta. The detection of Si XIII and Mg XI X-ray lines (Long et al. 2003) in 
the NW region of the remnant agrees with the hypothesis that the emission is coming from 
ejecta near the shock front. 

To find a proton temperature we have thus far used the width of the Ha line. However, 
the proton thermal speed is not simply equal to the velocity derived from the width of 
the Ha line at high temperatures. The cross section for neutral-proton charge transfer, 
the process that produces the broad Ha, falls off at high energies allowing for the neutral 
hydrogen distribution function to be narrower than that of the protons(Chevalier, Kirshner, 
& Raymond 1980). This results in an Ha profile that would incorrectly indicate a lower 
temperature than the actual proton temperature. 



4.3. Heating at the Shock Front 

From the current observations there are two ways to calculate the temperature of the 
ions. The first method to calculate the temperature is based on the thermalization of the 
bulk velocity of the shock. The second method uses the FWHM of the gaussian line fits 
as the thermal velocity that can be used to find the temperature. To determine the kinetic 
temperature of a species from the bulk thermalization we use the equation, 

3 

kTi = Y^ruiv^^^^^ (1) 
where the subscript i indicates the species, k is the Boltzman constant, T is temperature, 
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is the mass of the species and i^shock is the shock speed = 2890 km s'"^(GWRL02). This gives 
a temperature for O VI of 2.9 x 10^ K and for the protons of 1.8 x 10® K. The ratio of the 
temperatures is mass proportional, Toxygen=16Tproton, which is expected using this method. 
This heating occurs when some fraction of the energy of the shock speed is transferred to 
the thermal velocity of the protons or ions. 

Using the width of the O VI lines to calculate temperature, the temperature is measured 
to be 1.5 X 10^ K. The O VI temperature derived from the observed line width is less 
than that predicted by the kinetic temperature equation for no equilibration among particle 
species. The ratio of the temperatures indicates that O VI is heated to a temperature 48% 
less than the value predicted for mass proportional heating. Ions are being heated by a 
process other than the bulk fluid velocity thermalization or there is a heat loss mechanism 
for the ions. 

Heating of ions in coUisionless shocks has been studied by Berdichevsky et al. (1997) 
using heliospheric shock data. In examming O VII, it was found that the oxygen was pref- 
erentially heated 19-39 times more than the protons. In studying the solar wind, Lee & Wu 
(2000) assume greater than mass-proportional heating as part of the coronal heating process. 
As a consequence, ions non-adiabatically expand upstream (not being reflected by the shock 
front) and move with a velocity equal to their gyration velocity as they go upstream. These 
hot highly energized ions could act as a precursor that takes away a significant amount of 
energy. 

The current supernova observation of less than mass proportional heating lies in stark 

contrast to the heliospheric coUisionless shocks. Several factors and processes determine the 
extent of ion heating. The first comparison to be made is the speed of the shock relative to 
the local Alfvenic speed. The solar shocks propagate at 400-1000 km s^^. SN1006's shock 
is propagating at almost 3000 km s~^. The Alfvenic Mach number, the ratio of the shock 
speed to the square root average of the thermal and local Alfvenic speed, is < 10 for solar 
shocks but upwards of 200 for the supernova shock. The orientation of the magnetic field 
with respect to the normal of the shock front is also of importance as quasi-perpendicular 
shocks and quasi-parallel shocks are quite different. If the current magnetic field orientation 
for SN1006 is correct, the NW is propagating parallel to the ambient magnetic field, while 
both parallel and perpendicular shocks are observed in the solar wind. 

A measure of the importance of the magnetic field is the parameter (5. The plasma /3, 
the ratio of thermal to magnetic pressure upstream, is small (< 1) for hehospheric shocks. 

Using the parameters for SN1006 and the general value for the Galactic plane ISM magnetic 
field (~3 /uG), the NE has /? = 0.02 and in the NW {3 = 0.1. This indicates that the magnetic 
field pressure dominates over the thermal pressure at the ISM/remnant boundary as in the 
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solar wind, in contrast to the ISM which is assumed to have a /3 of unity. It is possible that 
the change in density from pre-shock to post-shock conditions is an important characteristic 
in the propagation and heating of ions in the coUisionless shock fronts. 

In order to determine the cause of the different heating found in the hehosphere and 
the supernova further investigation of the influence of pressure, density, Mach number and 
velocity on ion heating by the shock is necessary. 



4.4. Neutrals at the Shock Front 

In the analysis scheme used here from Chevalier & Raymond (1978) and Chevalier, 
Kirshner, & Raymond (1980), neutrals play a vital role. Neutrals undergo charge exchange 
or emit line radiation to produce the Ha and Ly-/3 emission. There are fast neutrals produced 
by the shock, as evident by the broad components of the Ha and Ly-/9 lines. This could 
create a neutral precursor for the shock(Smith ct al. 1994; Lim & Raga 1996). The hydrogen 
and oxygen neutral fractions are tightly coupled by charge transfer, thus information about 
the neutral fraction of hydrogen can be used to diagnose the neutral fraction of oxygen. 
These neutrals should become pickup ions like those seen in the solar wind (Vasyliunas & 
Siscoe 1976) when they pass through the shock and become ionized. Pickup ions like those 
in the hehosphere can then act as a high energy seed population for Fermi acceleration just as 
heliospheric pickup ions are the seed population for anomalous cosmic rays(Fisk, Kozlovsky 
& Ramaty 1974). 

The He II 4686A hue can be used as an indicator of neutral fraction due to its in- 
sensitivity to pre-shock neutral fraction and electron-ion pre-shock temperature equilibrium 
(GWRL02). In the NW, observations have shown He II emission lines (Raymond, Blair, Long 
1995)(GWRL02). The ratio of He I/He II can then be used to find an H neutral fraction 
which is a parameter in the relation of the Ha two component intensity ratio, Ibroad/Inarrow, 
and the electron-ion temperature ratio. In addition, GWRL02 calculated the pre-shock H 
population to be 90% ionized but the pre-shock He population is 70% neutral. Using the 
Ha broad-to-narrow intensity ratio calculated for 90% pre-ionized medium, the temperature 
ratio, Teiectron/Tproton, was fouud to be < 0.07, showing little to no equilibration between 
protons and electrons. Using the ratio of Teiectron/Tproton, we find an electron temperature 
of < 1.2 X 10^ K, approximately 1 keV, which is an upper limit that agrees with the value 
found by Long et al. (2003) of Tdcctron < 0.6 keV, but significantly less than the oxygen and 
proton temperatures found for this observation (1.5 x 10^ K and 1.8 x 10^ K, respectively). 
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5. Summciry 

In summary, the two shock regions of SN1006 studied here provide a unique cosmic 
laboratory for shocks and their acceleration processes. Clearly, the properties of the inter- 
stellar medium play a crucial role in shaping these shocks. We conclude with the following 
summary of our observations and interpretations. 

1. The material that the NE shock front is encountering is less dense than in the NW 
region, with a ratio n(NW)/n(NE)>4, and is best seen in the X-ray or radio wave- 
lengths. The NW shock front could be moving into a diffuse H I cloud or similarly 
dense region. 

2. The O VI line width of the NW shock indicates that oxygen ions are heated to a 
temperature less than 48% of the value predicted by mass proportional heating. This 
differs from the observations of other non-radiative coUisionless shock fronts such as 
those in the heliosphere which found ion temperatures 20-40 times in excess of the 
values predicted by mass proportional heating. The roles of density, pressure, magnetic 
field orientation with respect to the shock normal, velocity and Mach number should 
be examined to better determine the ion heating mechanisms. 

3. The plasma at the shock front has not had time to come to equilibrium via Coulomb 
collisions. The plasma is in a non-equilibrium state with energy distributed differently 

between species of the plasma. This is in agreement with the work on temperature 
equilibrium done by Ghavamian et al. (2002) who found the ratio of proton to electron 
temperature to be < 0.07 indicating a plasma far from equilibrium. 

The role of the neutral fraction of the ISM population in the charge exchange interaction 
should be examined in detail as it may greatly affect the outcome of the shock-ISM interac- 
tion. The rate at which the plasma becomes isotropic, the particle distribution, and the time 
scale to reach isotropic and Maxwellian conditions are in need of examination to understand 
the heating process present in the coUisionless shock. Further work will be done to model 
the plasma conditions in coUisionless shock fronts to include neutral fraction as well as ex- 
amine the role of electron population on O VI formation. This work should help advance the 
understanding of the shock acceleration of particles and the physics of a coUisionless shocks 
in varying environments. 

This work is based on observations made with the NASA-CNES-CSA Far Ultraviolet 
Spectroscopic Explorer. FUSE is operated for NASA by the Johns Hopkins University 
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Fig. 1. — Ha image of SN1006 taken by the CTIO Schmidt telescope from Winkler et al. 
(2003). Closeup images of the observed filaments are shown in the insets. In each case the 
interior box, drawn to scale, shows the location of the 30" x 30" FUSE LWRS aperture. The 
NW blowup is from the same Ho: image while the NE blowup is a 0.3-8 keV Chandra image 
(Long et al. 2003). 

Fig. 2.— Raw FUSE spectra from the Northeast and Northwest region of SN1006. The 

NW is offset from the NE by 0.5 for clarity. Geocoronal lines are marked. The Ly-/? peak 
dominates in both spectra. The NE is consistently fainter than the NW, but the geocoronal 
line intensities are similar. 

Fig. 3. — FUSE spectra from the NW, binned at 0.1 A with the NE subtracted as back- 
ground. The dotted dashed lines are the regions of the spectra that were excluded from the 
fits. The dashed lines are the fits for the Ly-/3, O VI 1032 and 1037 A lines, with FWHM of 
2290 and 2100 km s~^ respectively. The solid line represents the addition of the fits of the 
three spectral lines. 
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Table 1. Summary of UV Emission Lines in NW Filament of SN1006 



Ion 


Intensity 


Filament Length 


FWHM 


Temperature 


mion/mpT 


% 




(photons cm^^s^^arcsec"^) 


(arcsec) 


(km s"-') 


(Kelvin) 


(Kelvin) 


Mass 




(xlO-4) 




Observed 


from FWHM 




Prop 


H-ai 


2.1 


51 


2290 ± 80 


1.8 X 10* 






Ly-13 


4.0 


30 


2290(fixcd) 








He II^ 


0.99 


197 


2558±618 


5.7 X 10** 


7.2 X 10** 


79% 


C IV^ 


1.7 


197 


2641± 355 


1.8x 10^ 


2.2x 10" 


82% 


O VI 


3.1 


30 


2100 ± 200 


1.5x 10" 


2.9x 10" 


52% 


O VII^ 




60 


1775±261 


1.1 X lO'' 


2.9x 10" 


38% 



''Temperature derived from shock speed of 2890 km s ^ . 

References. — (1) Ghavamian et al. 2002; (2) Raymond et al. 1995; (3) Vink et al. 2003 



